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Abstract 



^ ' In this paper we consider black p-branes on square torus. We find an indication of 

H \ 

\ a cascade of Gregory-Laflamme transitions between black p-brane and (p — l)-brane. 

Through AdS/CFT correspondence, these transitions are related to the breakdown 
of the ^7(1) symmetry in super Yang-Mills on square torus. We argue a relationship 
between the cascade and recent Monte-Carlo data. 
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1 Introduction 



In general relativity a black string becomes unstable and changes to a black hole when the 
circle it winds becomes large compared with its horizon. This is known as the Gregory- 
Laflamme transition (see for review H). Similar instability exists also for neutral 
black p-brane wrapped on torus. In the case of square torus, in which all compactification 
radii are the same, usually only a black hole and a non-uniform brane are considered as 
candidates of the final state, probably because they respect the symmetry of the back- 
ground spacetime. However, there is no reason for adhering them - isotropy can break 
spontaneously, and we should also consider a decay from p-brane to {p — l)-brane. 

A motivation for considering such a decay comes from AdS/CFT correspondence 0, 
I, g. Through AdS/CFT, the Gregory-Laflamme transition is related 0, g |, [ig, [IT], || 
to the condensation of the spatial Wilson loop winding on compact direction Q in dual 
super Yang-Mills theory. Recently, Narayanan, Neuberger and collaborators have studied 
phase structure of pure bosonic Yang-Mills theory on torus intensively [T^|. According 
to their result, the spatial Wilson loops condensate not simultaneously but one- by-one 
even if compactification radii are the same. We can expect similar pattern of successive 
phase transitions in bosonic Yang-Mills on T" with 10 — n adjoint scalars, which is the 
high-temperature limit of the super Yang-Mills on T""*"^. Then, it is plausible that dual 
supergravity theory also has a cascade of Gregory-Lafiamme transitions. In this paper, 
we compare the free energies for several solutions in (super)gravity using a certain ap- 
proximation. The result suggests that such a phase structure exists. In IIA supergravity 
the free energy of Dp-braneQ becomes smaller than that of D(p— l)-brane at some critical 
temperature tc{p), which satisfies tc{p) > tc{p-i)- In the case of the transition from Dl- to 
DO-brane, the Gregory-Laflamme temperature appears slightly below tc'(i). This implies 
that the Dl-brane collapses to the DO-brane. We expect a similar pattern for any p. In 
Einstein gravity, we confirm this pattern explicitly for the neutral black branes (see Table 
|, Tableland Figure |l|). 

The organization of this paper is as follows. In §0 we consider the neutral black branes. 
Approximating compact directions transverse to brane with noncompact ones, we find an 
indication that p-brane does not decay to 0-brane directly, but its dimension decreases 
one-by-one, that is, p-brane decays to {p — l)-brane, then to {p — 2)-brane, and so on. 
In §1^, we summarize basic properties of weakly-coupled Yang-Mills theory, and then we 
study the phase structure of dual supergravity in We find similar phase structures in 
both sides. This result might be an evidence of AdS/CFT, although the coupling region 
we considered are different in both sides. 55 is devoted to discussions. 



■^In the following we call it simply "spatial Wilson loop" . 
"^Precisely speaking, smeared DO-branes obtained via T-dual. 



1 



2 Neutral black branes on M^"^ x 

In this section, we calculate the free energy of neutral black p-branes winding around 
square torus T" and evaluate critical temperatures tc{p), above which {p — l)-brane has 
smaller free energy than p-brane. We also calculate the Gregory-Laflamme (GL) tempera- 
ture tGL{p)i above which p-brane becomes unstable. In these calculations, we approximate 
compact directions transverse to brane by noncompact ones because of the lack of analytic 
solution for black branes in compact space. We find that 

■ ■ ■ < tc(p) < iGL{p) < ic(p-l) < ^GL(p-l) < • • • (2.1) 

for D < 12 and D — p > 4 which indicates the existence of a cascade of first order 
transitions^. 



2.1 Preliminary 

Let us consider the D-dimensional Einstein-Hilbert action in the asymptotically flat space, 

1 



d^x^/^R. (2.2) 



I EH 



We take the background spacetime to be R-^"*^ x T^, where compactification radii of T'^ 
are L for all directions. Suppose that a Schwarzschild black hole is placed in this space. 
If its horizon is much smaller than L, then we can approximate it using a Schwarzschild 
solution in M^. If the size of the horizon becomes comparable to L, the black hole cannot 
fit in and can wind on some cycles in T". This is the Gregory-Lafiamme transition. 
Preferred value of p can be determined by comparing the free energy. 

Let us begin with the simplest case, p = n. In this case, the metric can be written as 



where represents the horizon radius and < yi < L {i = 1 . . . n) . By requiring the 
regularity of the above metric, the Hawking temperature is given by 

r„) D — n — 3 / .N 

Th = r^. 2.4 



^In the case oi D > 13, (2.1) does not hold. For example, at = 13, we obtain tQi^^p-^ < t(j(^p-j for 
p < 3. Although this is consistent with |^9| , however, our approximation becomes subtle at such large 
dimensions as pointed out in . More careful treatment is needed in order to show whether or not our 
approximation could be applicable at larger dimensions. 
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The ADM energy (mass) and entropy of the black string ( |2.3| ) are given by [18 



M(n) = ^—j^{D -n- 2)no^^^2{rf)''~^~^L^, (2.5) 



1 



S{n) = -—^Qn-n-2{rP)''--'L\ (2.6) 



For p < n, because exact metric is not known, we approximate the transverse compact 
directions by noncompact ones. Then, the above expressions can be used by replacing n 
with p. In several examples, this approximation is known to be good even if the horizon 



is comparable to the compactification radius; see e.g. |26|. In the present case, the ratio 



vh/L is around 0.3 ~ 0.4 and we can expect that this approximation is valid. 



2.2 j9-brane vs. {p — l)-brane 

In order to determine the phase structure, let us compare the free energies for p-branes 
and {p — l)-branes fixing temperature Th and varying L. Then, because the Hawking 
temperature is related to the size of the horizon as 

'^H^ ~ Tl — ~ ■ iiidependent of p, (2.7) 



we have the relation 



The free energy of p-brane is given by 



When F{p) > F{p — 1), the {p — l)-brane is more favorable. Using (|2.9|) and (2^), this 
relation becomes 

nn r f D - p - 2\'^ Vt D -r,-l D - p - 

As shown in Table ^ the critical value tc{p) is the decreasing function of p. 

We can also compare the entropy with ADM energy fixed. Similar pattern of instability 
appears also in this case (see Table p. 

The argument above does not exclude the possibility that several phases co-exist . That 
the instability really arises can be confirmed by calculating the Gregory-Laflamme (GL) 
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mass fixed 


temperature fixed 


tc{i) 


1.27 


1.28 


tc(2) 


1.15 


1.17 




1.01 


1.04 



Table 1: The critical values of the circle for MJ x T^. 



critical temperature |I|], where ]?-branes become unstable thermodynamically |1^, |I5|. As 
we will show in the following, there is a relation^ 



(2.11] 



for D = 10. (Note that this relation does not hold for large D.) The relation ( |2.11|) 
suggests that the p-brane decays to the {p — l)-brane at t = tciip) as depicted in Figure [l|. 
Then it is plausible that the transition is of first order. 

Now let us evaluate tGL(p)- By approximating the compact directions transverse to p- 
brane by noncompact ones, the problem reduces to the evaluation of the GL temperature 
for p-brane wrapped on in R^-P X T^. This has been studied already in p9|, |23, ^ 
(see Table 0). 



^GL(l) 


1.30 


tGL{2) 


1.20 


tGL{3) 


1.08 



Table 2: The Gregory-Laflamme critical temperature for 



T3. 



Furthermore, we expect p-brane decays not to 0-brane directly, but to {p — l)-brane. 
We check this below for D = 10, p < 3 hj comparing free energies at tcLip)- Figure ^ 
shows that only [p — l)-brane has smaller free energy than p-brane at tciip)- This result 
makes us confirm the above statement is trueQ. Similar statement can be seen in ||30| , 
where n-brane on T" has been studied and it was shown that tachyonic mode develops 
along one of the n-directions 0. It is important to study instability around p-brane {p < n) 
in T" by metric perturbation and determine whether the tachyonic mode grows along one 



The authors arc grateful to J. Marsano for instructive comments concerning the derivation of ( 2.11 



^ Similar cascade can be found when we compare the entropies of p-branes with fixed mass. For 
example, for D = 10, n = 3, we have 5(0) < S{3) < S{1) < S{2) at t = tQ^^^y Therefore, 3-brane can 
decay only to 1- and 2-branes. Entropically, 2-branc is most favorable. Once 3-brane decays to 2-brane, 
then we have S{2) < S{0) < S{1) at t = tGi(2)i and hence 2-brane can decay to 0- and 1-branes. Again, 
1-brane is more favorable entropically. This result is consistent with the result in [ p9| . 

® We are grateful to E. Sorkin for informing us |3Q|. 
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Figure 1: A plausible picture of a phase transition. [Left] At t > tc{p); p-brane has 
smaller free energy and is favored. [Center] At t = tc(p), free energies of p- and {p — 1)- 
branes become the same. [Right] At t = tGL{p), p-brane becomes unstable and decays to 
{p — l)-brane. 



of the directions also in this case. Also it is necessary to calculate critical temperatures 
tc{i) and tcL{i) more precisely, because a few percent error can destroy the cascade [^ . 






Figure 2: Comparison of free energies at tc(p) in unit of ^/IGnG^ for D = 10, p < 3. 
Only [p — l)-brane has smaller free energy than p-brane and is favored. [Left] t = tc{3)- 
[Center] t = tc{2)- [Right] t = tc(i)- 

Alternatively, we can compare the entropy while taking the mass of p- and p — 1-branes 
to be equal. In this M decrease through radiation, S{p) becomes smaller than 

S{p - 1) at 



LD~3 (D-P- 2)(^-p-2)' (1]d„p„2)''"P~' 

IGttG^ {D-P- l){D-p-l){D-p-3) ^^^_^_^)D-p-3- 



(2.12) 



This M|p^p_i is monotonically increasing with p, and hence transitions can take place 
repeatedly. 
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3 Super Yang-Mills theory on T^+i and its phase struc- 
ture 

In this section, we consider the U{N) super Yang-Mills theory on Euclidean torus T""^^ (n = 
1, 2, 3).0 The action is given by 

Sym = ^ / dtl dx'---! dx" Tr {F^^ + 2{D,<p'Y - [<j)' , (p-"]^ + (fermion)} , 

(3.1) 

where ^ = 0, ■ ■ ■ ,n, x^ = t, (p^ {I = 1, ■ ■ ■ ,9 — n) are adjoint scalars, and fermions are 
anti-periodic along the thermal cycle. 

At high-temperature, KK modes along thermal cycle decouple and the model can be 
described by bosonic Yang-Mills on T". If L is small, spatial KK modes decouple and 
the model reduces to BFSS matrix model When both thermal and spatial cycles 

are small, it reduces to a bosonic version |^ of the IKKT matrix model [^. In the 
following, we determine the parametric region where the above reduction is justified, and 
then discuss phase structure there. 

By rewriting the action using dimensionless quantities, we have 

AT r^/tH fi fi 

SYM = -rT^ dt dx'--- rfx" Tr {F^^ + 2{D,<f)')^ - [0^ <p']' + (fermion)} , 

Jo Jo Jo 

(3.2) 

where 

A' = AL^-", tH = ThL, (3.3) 

and fields are also redefined to be dimensionless e.g. as A^^^ = LA°^'^. If spatial KK 
modes decouple before temporal KK modes do, then the action can be rewritten as 

SyM = Jyj dtTr{F'^, + ---} = - j rftrr|F2, + ...|, 

(3.4) 

where t = X'^^^t and F = X~''^^^F. Therefore, temporal zero modes decouple if 

tH > A'l/l (3 5) 

In the same way, if temporal KK modes decouple before spatial KK modes do, then 
spatial KK modes decouple when 

tH « ^. (3.6) 



The argument below is a generalization of that in [O] to n > 1. Here we take large- TV limit. 



6 



In summary, SYM on T^+"- can be described by bosonic IKKT model (i.e. can be 
reduced to zero- dimensional bosonic model) if 

A'^/=^ « tH « ^. (3.7) 

In this region, both Polyakov loop and the spatial Wilson loops have nonzero expectation 
values. At t// > ^, the model reduces to bosonic YM on with (10 —p) adjoint scalars, 
which can be studied through lattice simulation. In this region, the Polyakov loop has 
nonzero expectation value and hence the model is in the black hole/black string phase. 
Condensation of spatial Wilson loops takes place at t/^ ~ p-. 

Let us consider the case of n = 2.0 Then, at t// > ^, SYM is well approximated by 
two-dimensional YM with 8 adjoint scalars. Although a simulation of this model has not 
yet performed as far as we know, we can guess its properties from simulations on similar 
models. In |]16|, bosonic pure U{N) YM on and have been studied. By taking the 
sizes of one or two directions to be zero, they become bosonic YM on with one or 
two adjoint scalars, respectively. Let us take two compactification radii of to be the 
same value L. Then, the action is invariant under the exchange of two directions. Let the 
spatial Wilson loops to be 



Pexp i / dx^A 







(/i = l,2), (3.^ 



where the indices fi are not contracted in r.h.s. Spatial Wilson loops have nonzero expec- 
tation values if the global ?7(1)^ symmetry 

A^~^Af, + c^- In (3.9) 



is spontaneously broken. As discussed in |T^, the U{iy symmetry is actually broken and 
(W^) behaves as 

{Wi) = {W2) = {L>^L^}'^), (3.10) 
(1^1)^0, (1^2) = or (1^2)7^0, (1^1) = (^L(2) < L < ^L«),(3.11) 
(1^1)^0,(1^2)7^0 {L<^Lf^). (3.12) 

We can expect that this pattern does not depend on the number of adjoint scalars. We can 
also expect that the result is similar also for other values of n. Note that these transitions 
are likely to be of first order. 

Suppose that such a pattern of phase transition persists to strong coupling region. In 
the dual gravity theory, the eigenvalue distribution of spatial Wilson loops is related to 
that of DO-branes in T-dual picture, and condensation of spatial Wilson loop corresponds 



1 case has been studied in ||Tl| . Simulation results can be found in ||Tl|, |2^, |2^, |2j, |2£ 
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to the Gregory-Laflamme transition in DO-brane system ^ |10|, |TT], [T^. Therefore, 
transitions in super Yang-Mills theory can be regarded as a sequence of Gregory-Lafiamme 
transitions in gravity side. In the next section, we show such successive transitions do 
exist in supergravity. 



4 Cascade of Gregory-Laflamme transitions in IIA 
Supergravity 

In this section, we study the dual gravity theory of Dp-branes on T" corresponding to 
the super Yang-Mills theory on square torus discussed in the previous section and find an 
indication for the cascade. As we will show, the transitions take place in the parametric 
region where the T-dual picture, a system of DO-branes, is valid. 



4.1 Generalities 

We consider the near-extremal Dp-brane solution on T"^ which is dual to the super Yang- 
Mills theory on T" in finite temperature. As in §0, for p < n we approximate transverse 
compact dimensions by noncompact ones. 

In the near horizon limit the metric and dilaton are given by [l^ 



U 2 



9YM\/dpN 



tt7~P\ P 



i=l 



QYAI^/dpN 



3-p 



8-p (y 



2-V..2 f 9YMdpN\ 4 ^ ^^7_2p^2^^ f7 -p 



(4.1) 



(4.2) 



where < ?/j < L (i = 1, ■ ■ ■ , p) denotes the direction of T^. The Yang-Mills coupling 
constant is related to the string coupling at the spatial infinity Qs = e'^°° as 



(4.3) 



Requiring the regularity of the metric around the horizon th-, the Hawking temperature 
is determined to be 



H 



An^ydpX 



(4.4) 
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where A = Qym^ ■ Using the Einstein frame the ADM energy and the entropy of this 
solution can be computed as p!3| 



Sdp 



9-p 



2ii-2P7r^r(^)A2 



9-p 



28-2Pvr^r(^)A2 



(4.5) 
(4.6) 



We can use the metric ( |4.1| ) as a supergravity solution when the string excitation and 
winding modes wrapping on are much heavier than KK modes along S*^"^: 



\/AL 



yAL(5-p)/2 



1 « A {p = 3) 
<rjf<A^ (p = 0,l,2). 



(4.7) 
(4.8) 



T-dual picture 



By taking a T-dual in all the compact directions of T^, the metric and dilaton become 



ds =a' < 



U 2 



de + ^/d^u'-^dn'^ 



8-p 



+ 



^/dpX 

U 2 



dU 



J— + dy} 



f/7-p 



i=l 



= (2 



TT 



2 A / (ipA \ ^ 1 



(4.9) 



(4.10) 



where denotes the T-dualized circle coordinate and takes the rangeQO < yi < (27r)2/L. 
The metric after T-dual represents the uniform distribution of DO-branes on T^. The 
Hawking temperature, ADM energy and entropy are the same as those before taking 
T-dual, (H), §1) and 



Winding modes are negligible when 
ThL < 1, 

and the condition that the a' correction can be neglected does not change: 
Th<A^/(=^-p) (p = 0,1,2), KA (p = 3). 



We have used the relation between the dilaton before and after the T-dual e*^ = 
2ttR = L. 



(4.11) 



(4.12) 



-, where 
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4.2 Bp vs. D(j9 - 1) 



From ( |4.5| ) and ( [4.6| ), the free energy of Dp-branes is given by 

F(p) = J^-^ N^Ul^^'LP. (4.13) 

By using dimensionless parameters t = ThL, uq = UqL and A' = XL^~p, ADM energy, 
entropy and free energy can be rewritten as 



(9-p)5(p) 



p-3 2(7-p) 



E(p) = 2^'''' N'\^t^ , (4.14) 
Sip) = STiB{p)N^^'0t^\ (4.15) 
F{p) = -^^^^^N^Xpt^, (4.16) 



b — p 



2(7~P) 

5 — p 



Here (p) indicates that these quantities depend on p in general. Note, however, that the 
effective 't Hooft couphng A'^^^ is indeed independent of p. This can be seen as follows. 
We can write the 't Hooft coupling on using ( [4. 3D Q 



A(p) = ^?^MA^=(27r)-V^'a'^. (4.17) 



Since the string coupling on the and its T-dual is related as 



.(0) ap) / Va' 



we have 

A(p) = A(o)L^ (4.19) 
Then it is clear the effective coupling does not depend on p: 

^(p) = A' = A(o)j 

Temperature fixed comparison 



X' =X' = A(o)Ll (4.20) 



We compare F{p) and F{p — 1) at the same Th and L (and hence t) in the strong coupling 
region A' ^ 1. 



Using ([4.16|) we can easily see F{p) > F{p — 1) is realized when the dimensionless 



temperature t becomes lower than the critical temperature tc{p) 

M A(v) = ((i^mi^ 



10 



^C(4) 


2.92/yA' 


tc{3) 


2.87/v/A' 


tc{2) 


2.67/v/A' 


tc(i) 


2A0/VX' 



Table 3: The critical effective temperatures 



As shown in Table ^ tc(p) is the increasing function of p. 

These critical points lie in the parametric region ( [4.11| ) where DO-brane picture is 



valid. In this region KK modes along T" can also be neglected. Note that, in this picture, 
smaller t corresponds to larger compactification radius L ~ 1/L. Therefore, the relation 

tc{i) < ■ ■ ■ < tc{p~i) < tc(p) (4.22) 

indicates that for large L low dimensional object is favored similarly to the Schwarzschild 
case. 

Above we have assumed that the non-uniform phase where DO-branes collapse in 
some circle directions can be approximated by those on fiat noncompact space. In several 
examples, this assumption is known to be valid when the horizon is smaller than the radius 
of the circle. In the present case, the horizon of the DO-branes on and the radius of 
the T-dualized torus before taking near-horizon limit are given as r'f^ = u^^a'/L and 
L = (27r)^a'/L. Let us consider the ratio at the critical temperature tc{p) 



aip)^"-^ =77^ r"X^^^""M • (4.23) 



L 



1 / An^/d^A^p) 



2 

5-p 



(27r)2 I 7-p 



After the transition Dp-branes become D{p — l)-branes and the radius of the horizon 
changes to 



r 



(p-i) 



An^/d^A{p) 



2 

6—p 



hip) = = I - V I .4 24) 

Above ratios take values around 0.3 ~ 0.4 and we can expect that the assumption is valid 
(see Table ^ and Figure We can also see that a{p — 1) < h{p) for p < 4. This relation 
indicates that, after a transition from p-brane to {p — l)-brane, the horizon of {p — 1)- 
brane is larger than the next critical radius ^\t=tc{p-i) ^'^^ hence the next transition 
can take place when the horizon shrinks to ^^k=<c(p-i)- ■'■^ ^^^^ ^^Yi transitions take 
place repeatedly. 
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We return to the metric before taking near- horizon hmit. 
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p 


a{p) 


b{p) 


4 


0.302 


0.329 


3 


0.323 


0.369 


2 


0.351 


0.400 


1 


0.379 


0.427 



Table 4: The ratio between the horizon and the torus radius at the critical temperature. 



In the present case, we have not calculated the GL temperature tGL(p)- In ^gl(i) 
has been calculated, and the result is 



2.29 

torn - 

This is slightly lower than tc{i) as expected. If other tGL(p)'s satisfy 

^C(p-l) < icLip) < tc{p), 



(4.25) 



(4.26) 



then a cascade can take place as first order transitions similarly to the case of dual SYM. 
Such a cascade resembles also to the Schwarzschild case. 




Figure 3: The plot of a{p — 1) (red) and b{p) (green). 
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5 Discussion 



There are several directions for future studies. First of all, it is important to justify the 
approximation in this paper in which compact directions are replaced with noncompact 
ones, because phase structure is sensitive to the small changes of physical quantities; as 
discussed in [^|, a few percent error in the entropy can change the magnitude relation 
of critical temperatures and lead us to the wrong conclusion. At least in the case of D- 
brane discussed in §4, AdS/CFT correspondence suggests that this approximation gives 
the correct phase diagram. Secondly similar calculations for other compactifications are 
important. If compatified dimensions are curved, our approximation would not be valid at 
all, and more careful treatment would be necessary. In addition, it is necessary to consider 
other solutions such as a non-uniform string solution; at present, we showed that uniform 
p-brane solution is unstable and {p — l)-brane is (meta-)stable at a critical temperature, 
but it might be possible that p-brane decays to another stable solution. Also, it will be 
nice if it can be determined whether ( 4.26| ) holds or not. 

It is possible to perform a similar study in gravity theories with other boundary con- 
ditions. From the point of view of the super Yang-Mills theory, such a study will be useful 
to clarify the phase structure of SYM on Tp+^ [Hl-H 

In large- gauge theories, if the global f/(l) symmetry is not broken then expectation 
values of Wilson loops are independent of compactification radii through Eguchi-Kawai 
equivalence Using this property, simulation cost could be reduced to a large extent 
||T6| . Determination of the f/(l)-unbroken region is important also for this reason. 
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